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Abstract 



The energy dependence of chirally expanded irN isoscalar and isovector 
amplitudes bo(E) and b\ (E) respectively, for zero-momentum off shell pi- 
ous near threshold, is used to impose the minimal substitution requirement 
E — > E — V c on the properly constructed pion optical potential within a 
large-scale fit to 100 pionic-atom data across the periodic table which also 
include the recently established 'deeply bound' pionic atoms of Pb and Sn. 
This fit cannot be reconciled with the well known free-space values of the ttN 
threshold amplitudes. In contrast, introducing the empirically known energy 
dependence for on-shell pions leads to a better fit and to satisfactory values for 
the ttN threshold amplitudes. The difference between these two approaches 
is briefly discussed. 
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I. INTRODUCTION AND METHODOLOGY 



The recent observation of Is and 2p 'deeply bound' ir~ atomic states in isotopes of Pb 
[1-3] and very recently also of such Is states in isotopes of Sn [4] has triggered renewed 
interest in the issue of partial restoration of chiral symmetry in dense nuclear matter [5-13]. 
In a nutshell, it was argued that since (i) the pion in deeply bound states with relatively large 
neutron excess charts a fairly dense portion of the nuclear medium, and since (ii) the most 
influential term of the optical potential V opt for this situation is generated by the s-wave 
isovector nN threshold amplitude b±, and since (iii) b± in free-space is well approximated in 
lowest chiral-expansion order by the Tomozawa- Weinberg expression [14] 

&i = -£| = -0.08 m; 1 , (1) 

then deeply-bound states could yield valuable information on the dependence of f n on the 
density p. The pion decay constant /„- serves as an order parameter for the spontaneously 
broken chiral symmetry in hadronic physics, and its free-space value /„-= 92.4 MeV should 
go to zero in dense matter if and when chiral symmetry is restored. Indeed, it has been 
known for quite some time that the renormalized value of b\ required to fit pionic-atom data 
is about —0.12 m^ 1 [15,16] clearly more repulsive than the free-space value —0.09 m~ x [17]. 

Our most extensive recent work [11] has shown, however, that the deeply bound states 
by themselves on statistical grounds are insufficient to draw firm conclusions about whether 
or not bi is renormalized in dense matter. In fact, contrary to the expectation (i) above, 
the pion in deeply bound Is states does not chart higher- density regions of the nucleus 
than it does so in 'normal' Is states in light nuclei. It was shown in Ref. [11] that only 
by using a substantially larger data base that includes plenty of normal pionic atom data, 
and carefully considering uncertainties in the knowledge of neutron density distributions, it 
becomes possible to make a meaningful statement on the renormalization of the isovector 
threshold amplitude b±, i.e. &i = —0.108 ± 0.007 m~ x . It is tempting to ascribe this value 
for bi, using Eq. (1), to a renormalization of /„. in the nuclear medium. 

Recently, Kolomeitsev, Kaiser and Weise [12] have suggested that pionic atom data could 
be reproduced using a pion optical potential underlain by chirally expanded nN amplitudes, 
retaining the energy dependence of the amplitudes bo(E) and b\(E) for zero- momentum 
(q=0) pions in nuclear matter in order to impose the minimal substitution requirement 
E — > E — V c , where V c is the Coulomb potential. This has the advantage of enabling one 
to use a systematic chiral expansion as an input [18], rather than singling out the leading 
order term Eq. (1) for b\. Kolomeitsev et al. applied this programme to study the shifts and 
widths of the Pb and Sn pionic deeply bound states and reported substantial improvement 
in reproducing these data with only minimal phenomenological input, mostly limited to the 
p-wave irN interaction to which allegedly these data are insensitive [12,13]. One could argue, 
however, that the deeply bound states do not offer sufficient variation over the energy range 
spanned by the bulk data on pionic atoms which include both 'shallow' states as well as 
'deep' states. In the present work we test whether or not the energy dependence of this 
chiral expansion provides a satisfactory description for the bulk of pionic atom data. 

The data used in the present work consist of 100 strong-interaction shifts and widths, 
stretching from 20 Ne to 238 U [11]. The Klein-Gordon equation solved for the pionic-atom 
eigen energies is given by [11,19]: 
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V 2 - 2fi(B + V c ) + (B + V c f - 11(E)} i> = (h = c = 1) (2) 

where p is the pion-nucleus reduced mass, B is the complex binding energy and V c is the 
finite-size Coulomb interaction of the pion with the nucleus, including vacuum-polarization 
terms. The pion-nuclear polarization operator H(E) is given by the standard Ericson-Ericson 
form [20] 

IT = 2pK P t = q(r) + V • a(r) V , (3) 

with the s-wave part of V opt 

q(r) = -4tt(1 + -^){&o(r)[p n (r) + p p (r)] + h\p n (r) - Pp (r)}} 

-4tt(1 + ^)4 J B oPn (r) Pp (r) . (4) 

In these expressions p n and p p are the neutron and proton density distributions normalized 
to the number of neutrons N and number of protons Z, respectively, and M is the mass of 
the nucleon; q(r) is referred to as the s-wave potential term and a(r) is referred to as the 
p-wave potential term. The function bo(r) i n Eq. (4) is given in terms of the local Fermi 
momentum fcp(r) corresponding to the isoscalar nucleon density distribution: 

b~o(r) = bo-^(bl + 2bl)k F (r), (5) 

where the quadratic terms in 6 and b\ represent double-scattering modifications of bo. In 
particular, the b\ term represents a sizable correction to the nearly vanishing linear bo term. 
Similar double-scattering modifications of b±, as well as other correction terms to 11(E) 
listed in Refs. [12,18], were found by us to yield negligibly small effects and will not be 
further discussed below. The complex parameter B is due to s-wave absorption on pairs 
of nucleons. Its microscopic evaluation is outside the scope of chiral perturbation theory. 
Finally, the p-wave term a(r) is a standard one with the same form as in Ref. [11]. 

The chiral expansion of the ttN amplitudes for q = at the two-loop level is well 
approximated by the following expressions [12,18]: 



where a is the ttN sigma term, a ~ 50 MeV [21], is the nucleon axial- vector coupling 
constant, g\ = 1.27, f3 and 7 are tuned to reproduce the threshold values bo(m n ) ~ and 
bi(m n ) = — 0.0885lo;oo2i [17] respectively. For bo, in view of the accidental cancellations 
that lead to its near vanishing we limit our discussion to the f~ 2 term in Eq. (6), therefore 
choosing f3 = am' 2 . The next, f~ 4 term is much bigger than the scale of variation from zero 
expected for the threshold value and its inclusion here would appear somewhat dubious; if 
included, it would increase the energy dependence from the conservative estimate adopted 
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by us. Implementing the minimal substitution requirement in the calculation of pionic atom 
observables, the constant parameters 60,1 of the conventionally energy independent optical 
potential have been replaced in our calculation by 

b ,i(r) = b 0tl -6 0tl (RBB + V c (r)), (8) 

where 5 ,i — db 0> i(E) / dE is the appropriate slope parameter at threshold, ReB is the (real) 
binding energy of the corresponding pionic atom state and V c (r) is the Coulomb potential. 
The constant fit parameters 60,1 are then expected to agree with the corresponding free nN 
threshold amplitudes if the energy dependence is indeed responsible for the renormalized 
values found in conventional analyses. The added piece proportional to S in Eq. (8) is 
dominated by the attractive V c (r). Since the slope parameters S from Eqs. (6,7) are negative, 
this added piece is always repulsive, in agreement with Refs. [12,22]. 

Before testing the above 'chiral' energy dependence for off- shell q = pions we present 
results for the empirically known on-shell nN amplitudes, when the pion energy E and its 
three-momentum q are related by E 2 = m 2 + q 2 . This choice corresponds to the original 
suggestion by Ericson and Tauscher [22] to consider the effect of energy dependence in pionic 
atoms. Ericson subsequently [23] pointed out that, for strongly repulsive short-range NN 
correlations, the on-shell requirement follows naturally from the Agassi-Gal theorem [24] for 
scattering off non-overlapping nucleons. The corresponding nN amplitudes will be denoted 
below as 'empirical'. Figure 1 shows the energy dependence of the empirical bo(E) and 
bi (E) on-shell amplitudes as derived from the SAID data base [25]. The value of b (E) 
at threshold is very close to zero and the empirical slope So, which corresponds to adding 
repulsion in Eq. (8), is quite well-determined over the whole relevant energy range, having 
changed little since the classical KH80 analysis [26] of the pre pion-factories data to the most 
recent analyses of modern data. We note that the slope of the q = chiral bo(E) amplitude 
of Eq. (6) is larger than the on-shell empirical slope by about 60%. For the empirical b\{E), 
its value at threshold has also changed little since KH80 to the present day analysis, and 
the slope of the empirical amplitude is essentially zero, in contrast to the fairly large slope 
of the q = chiral amplitude of Eq. (7). 



II. RESULTS 

The present analysis is based on the 'global 3' data set of Ref. [11] consisting of 100 
data points from 20 Ne to 238 U. For the nuclear density distributions p p and p n we adopt 
the procedure of Ref. [11] where p p is obtained from the experimental charge distribution 
by unfolding the finite size of the charge of the proton, and where simple but physical 
parameterizations are used for p n . A key quantity in this context is the difference r n — r p 
between the root-mean-square radii. Relativistic mean field (RMF) calculations [27] yield 
to a good approximation [11] 

N ~ z 

r n -r p = a— - — + rj , (9) 

with the values a=1.51±0.07 fm, i]=— 0.03±0.01 fm. A similar expression, but with a = 1.0, 
was obtained by analyzing strong interaction effects in antiprotonic atoms [28]. Owing to 
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the strong correlation between the values assumed for r n — r p and the values of b\ derived 
from x 2 fits to pionic atom data, we have varied the neutron-excess parameter a over a 
wide range, with the expectation that a value in the range of 1.0 to 1.5 will represent on the 
average the 41 nuclei in the present data base. 

Figure 2 shows results for the 'conventional' model (left) for which 5 ,i = and for the 
'empirical' model (right) as function of the neutron-excess parameter a in Eq. (9) for two 
shapes of neutron densities. The dependence of the quality of fits on the shape of the neutron 
density distribution and the vanishingly small sensitivity of the derived values of b\ to this 
shape are demonstrated by using either the 'skin' or the 'halo' shape, as discussed in Ref. 
[11]. The results for the conventional (energy independent) model are practically the same as 
in Ref. [11] in spite of adopting now the 'current' SAID values [25] for the p-wave parameters 
Co = 0.21 m~ 3 and c\ = 0.165 m~ 3 , instead of the values 0.22 and 0.18 m~ 3 respectively 
used beforehand. This slight change was made for consistency, since the slope parameters 
8 for the empirical model were taken from the 'current' SAID analysis. In fact, we also 
incorporated in the right hand side (rhs) of the figure the SAID weak energy dependence of 
Co (ci is essentially energy independent). The results on the rhs of the figure show that, for 
the 'skin' shape of p n and with the introduction of the empirical energy dependence of the 
amplitudes, the minimum in the x 2 curve has shifted slightly towards the acceptable region 
of a =1.0 to 1.5, and the value of b\ for that minimum is in agreement with the free nN 
value marked by 'exp.' in the upper panels of the figure. It is self-evident that the 'halo' 
shape for the neutron density distributions cannot be reconciled with the data. Finally, we 
add that the resulting values of bo for both the conventional and empirical models are close 
to zero, well within the experimental error [17]. 

A comment on the 'anomalous s-wave repulsion' in pionic atoms is here in order. The net 
effect of the nearly vanishing parameter b , of the repulsive b\ and of the phenomenological 
parameter Rel?o has been known [19] to produce a repulsive potential inside nuclei which is 
twice as large as expected. This is due to the combined action of the too repulsive b\ and of 
ReBo which turns out to be too repulsive compared to the expectations that |Rei?o| <Imi?o 
(see also Ref. [7]). For the fits mentioned above we obtain for the conventional potential 
Im J B =0.053±0.002 m" 4 and Re^o = -0.10 ± 0.03 m~ A . Although the latter is determined 
to a moderate accuracy, we note that setting its value to zero while repeating the fits leads 
to a significant increase in the resulting x 2 value and to a value for bo which is incompatible 
with experiment (cf. Table 4 of Ref. [11]). Using the empirical bo,i(E) we find KeBo = —0.07 
± 0.03 m~ A . We conclude that using the empirical energy dependence the anomaly in Re_B 
is reduced, whereas there is essentially no anomaly in the parameter b±. 

Figure 3 shows results for the 'chiral' model, when either 6 (left) or b\ (right) is made 
energy dependent according to Eqs. (6) and (7), respectively. The left hand side of the figure 
shows that the quality of the best fit, upon incorporating only the energy dependence of 
the chiral bo amplitude, is significantly inferior to the corresponding best fit obtained using 
the conventional, energy independent model (shown in Fig. 2). Furthermore, the value of 
a at the x 2 minimum is unacceptably large and the corresponding value of b\ is in sharp 
disagreement with the experimental free 7iN threshold value. The rhs of the figure shows 
good fits with almost acceptable values for a and for bi, upon incorporating only the energy 
dependence of the chiral b\ amplitude, again for the 'skin' shape of the neutron density. 
However, incorporating the energy dependence of the chiral b (E) (even within the limited 
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scope of using only the f~ 2 term on the rhs of Eq. (6)) on top of that for bi(E), leads to 
substantial disagreement between the resulting best-fit value for b\ and the threshold value 
bi(m w ) which is marked by 'exp.' in Fig. 3. We conclude that, at present, the energy 
dependence generated by the chirally expanded s-wave nN amplitudes of Eqs. (6,7) fails 
badly in reproducing consistently the bulk of pionic atom data. 

III. DISCUSSION AND CONCLUSIONS 

In the present work we have demonstrated that the consistency between pionic-atom 
data and the free nN threshold amplitudes is greatly improved by using just the on-shell 
energy dependence of the nN s-wave amplitudes, in accordance with the original suggestion 
made by Ericson and Tauscher [22] . The idea behind using this empirical energy dependence 
is the same one as used for constructing the multiple-scattering series for short-ranged nN 
interactions occurring within an assembly of largely non-overlapping nucleons [20,23]. Mul- 
tiple scattering is naturally described in this idealized limit as occurring on shell. Whereas 
using off-shell q=0 pions in chiral expansions is motivated by the ground-state wavefunction 
description of pions in nuclear matter, applying this limitation to the construction of the 
pion-nuclear optical potential that generates pionic-atom wavefunctions is questionable. 

We have also shown that the energy-dependent chiral amplitudes given by Eqs. (6,7) for 
q=0 off-shell pions do not produce consistent or good global fits to pionic atom data. This 
conclusion is not at odds with the observation made by Kolomeitsev et al. [12] that the q=0 
off-shell chiral amplitudes work well, and with no need for a dispersive term ReB Q for the 
few deeply bound states available at present, since partial data sets of this kind do not have 
sufficient statistical significance to decide one way or another on this issue [10,11]. In fact, 
as good average reproduction of these deeply-bound data is reached within a wide classs 
of optical potentials, including our 'empirical' energy-dependent potential of the present 
study. We defer this and other ramifications of the present analysis for a forthcoming 
detailed publication. Given the fact that the on-shell nN amplitudes provide by far a better 
description of pionic atom data than the extremely off-shell q=0 chiral amplitudes do, we 
conclude that chiral dynamics is not yet at a stage of being tested in pionic atoms. 

Finally, it should be emphasized that we have strictly adhered in the present calculation 
to imposing minimal substitution, E — > E — V c , on the pion-nuclear polarization operator 
H(E) within the Klein-Gordon equation (2). Nowhere have we renormalized the threshold 
value of the nN isovector amplitude bi of Eq. (1) by renormalizing the pion decay constant 
f n — > f n (p) in dense matter [5]. This latter prescription which appears to be rooted in 
the underlying chiral symmetry has been discussed extensively in the context of pionic 
atoms [6-9,12,13], but according to Refs. [12,13] it need not be applied once the full energy 
dependence of II (E) is incorporated. 

We wish to thank Wolfram Weise for making critical remarks on a preliminary version 
of this manuscript. This research was partially supported by the Israel Science Foundation 
grant No. 131/01. 
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FIG. 1. ttN empirical s-wave scattering amplitudes as function of laboratory energy for on-shell 
pions from the SAID data base [25] . 
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FIG. 2. Pionic-atom fits using 'conventional', energy independent ttN amplitudes (left panels) 
and 'empirical', energy dependent SAID amplitudes (right panels) as function of the neutron-excess 
parameter a, Eq.(9), for two shapes of neutron densities. Lower part: values of x 2 for 100 data 
points from 20 Ne to 238 U. Upper part: best-fit values of b\ vs. the free ttN threshold value [17] 
marked by 'exp.' within the dotted horizontal lines. 
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